CONSTRUCTION OF DISCRETE ANALOGUE OF THE 
DIFFERENTIAL OPERATOR ^ + 2^ + 1 AND ITS 



(N| ■ PROPERTIES 

o 

(S| ■ A.R. HAYOTOV 

o 

Q 

in 



Abstract. In the present paper the discrete analogue of the differential 
operator -2-^ + 2-^ + 1 is constructed and its some properties are proved. 



The optimization problem of approximate integration formulas in the 
modern sense appears as the problem of finding the minimum of the norm 
of a error functional I given on some set of functions. 
^ , The minimization problem of the norm of the error functional by coeffi- 

' cients was reduced in [1] to the system of difference equations of Wiener- 

Hopf type in the space where is the space of functions with 

square integrable to— th generalized derivative. Existence and uniqueness of 
the solution of this system was proved by S.L. Sobolev. In the work [1] the 
description of some analytic algorithm for finding the coefficients of optimal 
. formulas is given. For this S.L. Sobolev defined and investigated the discrete 

' analogue of the differential operator of the polyharmonic operator A m . The 



problem of construction of the discrete operator D)J [p\ for n— dimensional 



^ case was very hard. In one dimensional case the discrete analogue D^[(3] 

of the differential operator ^ m was constructed by Z.Zh. Zhamalov and 
Kh.M. Shadimetov [2,3]. 

Further, in the work [4] the discrete analogue of the differential ^ x ™ m — 
fjim-i was constructed. The constructed discrete analogue of the operator 
y\ \ j x '2m — dx' 2m ~' 2 was a PP ne d for finding the coefficients of optimal quadrature 

formulas in the space 

Here, we mainly use a concept of functions of a discrete argument and 
the corresponding operations (see [1]). For completeness we give some of 
definitions. 

Let \f$\ = h/3, /3 G Z, h = jj, N = 1,2, ... . Assume that ip(x) and ip(x) 
are real-valued functions defined on the real line M. 

Definition 1. The function (p[(3] is a function of discrete argument if it 
is given on some set of integer values of f3. 
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Definition 2. The inner product of two discrete argument functions 
ip[/3] and tp[P] is given by 

oo 
v=—oo 

if the series on the right hand side converges absolutely. 

Definition 3. The convolution of two functions </?[/?] and ip[j3] is the 
inner product 

oo 

M*^\fi] = [v>bU\fi-i\]= E vb]-i>[P-i\- 

7=— oo 

In the present paper we consider the problem of construction of the dis- 
crete function D[(3] which satisfies the following equation 



D\p]*1>\J3]=6\p], (1) 



where 



G'[/3] = ^^(sin[/3]-[/3].cos[/3]), (2) 

5[/3] = | q' _^ q' is the discrete delta function. 

The discrete function D[/3] has important role in calculation of the coef- 
ficients of the optimal quadrature formulas in the Hilbert space 

^2(^2) = jv 9 : [0, 1] — > M <p' is absolutely continuous and ip" € £2(0, 1) j, 

equipped with the norm 



\<p\\ = { J (P2 (d/dx) ip(x)) 2 dx 
.0 

where P2(d/dx) = d 2 /dx 2 + 1. 

The equation (1) is the discrete analogue of the following equation 

(^ +2 ^ + 1 ) G(x) = S(l) ' (3) 

where G(x) = S1 §°( x ) ( s i n x —x- cos x) , 8(x) is Dirac's delta function. Moreover 
the discrete function D[f3\ has similar properties as the differential operator 
lx% + ^Ix 1 + ^' *h e zeros °f the discrete operator D[/3] are the discrete 
functions corresponding to the zeros of the operator ^ + 2^ + 1. The dis- 
crete function D[(3] is called the discrete analogue of the differential operator 

dx 4 ^ dx 2 ' 

The main results of the present work are the following theorems 
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Theorem 1. The discrete analogue of the differential operator Jpr + 
+ 1 satisfying equation (1) has the form 



m = - 2 



Ax-Af 1 " 1 , |/3|>2, 

l + ^i, , 101 = 1, (4) 

2h cos 2fo— sin 2h , ^4j_ o ri 

sin /1— A cos h "T" Ai ' ^ — u ' 



sin h — h ■ cos /i 
w/iere 

^ _ 4/i 2 sin 4 /i • A? _ 2/i - sin 2h - 2 sin h ■ \J K 1 - sin 2 h 

1 (A 2 - l)(sin/i - /icos/i) 2 ' 1 2(/icos /i - sink) 

I Ai| < 1, h is a small parameter. 

Theorem 2. The discrete analogue D[/3] of the differential operator Jjt + 

2^j + 1 satisfies the following equalities 
1) D[/3] * sm[/3] = 0, 
£>[/?]* cos [/3] = 0, 
s;D[0]*[0]sin[/3] = O, 
^ £>[£] * \j3] cos[/3] = 0, 

5) D[(3]*m = m- 

Here G[/3] is defined by (2), and 5[(3] is discrete delta function. 

In the proofs of these theorems we need the following well known for- 
mulas from the theory of generalized (distribution) functions and Fourier 
transformations (see, for instance, [1]) 

00 00 

F[<p(x)] = J V {x)^dx, F-^p)} = J tptfe-^dp, (5) 
—00 —00 

F[ip *ip] = F[cp] ■ F[ip], (6) 

F[ ( p-^]=F[ ( p]*F[^], (7) 

F[S^(x)] = (-2mp) a , F[S(x)] = l, (8) 

S(hx) = h-^ix), (9) 

5{x-a)-f(x) = S(x-a)-f(a), (10) 

6^\x)*f{x) = f^\x), (11) 

00 

Mz)= T, 8(x-P), ^e 2 ™^ = ^5(x-/3). (12) 

0=-oo P 13 

Proof of Theorem 1. 

According to the theory of periodic generalized functions and Fourier 
transformations instead of the function D[/3] it is convenient to search the 
harrow-shaped function (see [1]) 

00 

D(x)= Yj D[p]5(x - hp). 

/3=-oo 
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In the class of harrow-shaped functions equation (1) will be in the follow- 
ing form 

D Or)* G (x) = 5(x), (13) 

1-7 00 

where G (x) = G[(3]5(x — hj3) is the harrow-shaped function corre- 

/3=-oo 

sponding to the function G[jS\. 

Applying the Fourier transformation to both sides of equation (13) and 
taking into account (6), (8) we have 

F[D {x)\ = 1/F[G (x)]. (14) 

Now we calculate the Fourier transformation F[G (x)]. Using equalities 
(10) and (12), we get 

G {x) = h~ l G{x) ■ ^{h^x). (15) 
Further, use of equation (9) gives 

FiMh^x)} = hfo{hp). (16) 
Then, taking account of (15), (16) and (7), we obtain 

F[G (x)] = F[G(x)]*Mhp). (17) 

To calculate the Fourier transformation F[G(x)] we use equation (3). Taking 
into account (11), we rewrite equation (3) in the following form 

(5 (4) (x) + 25 {2) {x) + 5{x)) * G(x) = 5(x). 

Hence, keeping in mind (6), (8), we have 

F ^ = ^-D^+ir <18) 

Taking into account (18) from (17) we get 

f h 3 ^ 1 

Then from (14) we have 
(2vr) 4 



F[D](p) = 



£ (P-h{p + ±mp-h( P -±)y 



(19) 



/3=-oo yH yl "2n>> W 2n' 

Suppose the Fourier series of the function F[D](p) has the following form 

oo 

F[D](p)= D[p]e 2mph ^ (20) 
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where D[/3] is the Fourier coefficients of the function F[D](p), i.e. 

h- 1 

D[(3] = j F[D](p) e- 2mph "d P . (21) 
o 

Applying the inverse Fourier transformation to both sides of (20) we obtain 
the following harrow-shaped function 

oo 

D(x)= D[/3]5(x - hp). 

/3=-oo 

Then according to the definition of harrow-shaped functions the discrete 
function D\j3] is searching function of discrete argument D[f3]. Here for 
finding of the function D[f3] we will not use the formula (21). We will find 
it by the following way. 

To calculate the series (19) we use the following well known formula from 
the residual theory (see [6], p. 296) 

oo 

E /(/ 5 ) = - E res(7rctg(7r*) ■/(*)), (22) 

0=— OO Z\,Z1,...,Z n 

where z±, Z2, z n are poles of the function f(z). 

We denote f(z) = -. — , , , i — r? rrw- Here z\ = hip + and 

Z2 = h(p — ^) are the poles of order 2. Then taking into account (22) from 
(19) we have 

(Svr) 4 ' 



h 3 



^ res(7TCtg(7T2!) • f(z)) 



Zl,Z2 



F[D](p) 
Since 

2vr 4 / h sin(2vr/ip + h) 

res (vrctg(7rz) • f(z)) = — — - — — + 



(23) 



z=zi h 3 \ 1 — cos(2irhp + h) 1 — cos(2irhp + h) 

. . , ,. 27r 4 / /i sin(27r/ip — h) 

res (7rctg(7T2;) • f{z)) = —-rr- ; — ; rr ; — ; rr i ■ 

z=z 2 y BV ' Jy " h 3 \1 - cos(2nhp - h) 1 - cos(2nhp - h) J 

Denoting A = e 27Tiph , using the last two equalities and taking into account 
the following well known formulas 

e zi _|_ e ~ zi 
cosz = , sinz = 



2i 

after some simplifications from (23) for F[D](p) we get 

2 A 4 - 4 cos h A 3 + (2cos(2/i) + 4) A 2 - 4 cos h A + 1 



F[D]( P ) 



Smh-hcosh \ ( \2 , 2fe-sin(2fc) 



\ ( \2 i 2fe-sin(2fe) , 
A ^ A + sm/i-/icosh A + 1 



(24) 
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To find the explicit form of the discrete operator D[f3] the equality (24) we 
expand to the sum of elementary fractions. Since the polynomial Q2W = 

X2 + sLh-h { cSh X + 1 has two real roots 

2h - sin(2/i) - 2 sin h ■ \Jh 2 — sin 2 h 
Ai = 



A 2 



2(/icos h — sin/i) 
2h - sin(2/i) + 2 sin h ■ \J h 2 - sin 2 h 



2(/icos h — sin/i) 

and Ai ■ A2 = 1, |Ai| < 1. 
Then from (24) we have 

2 A 4 - 4 cos h A 3 + (2cos(2/i) +4)A 2 -4cos/i A + 1 



sin h — h cos h 



A + sin h- h cosh A + L ) 



2 / 2/icos(2/i) -sin(2/i) A Ai Si 
'A + . ' ^ + ^ + -, — + 



sin h — h cos h \ sin h — h cos /i A A — Ai A — A2 

'(25) 

For finding unknown coefficients A, A\ and B\ in the equation (25) we put 
A = 0, A = Ai and A = A2. 
For A = 

A = l, (26) 

for A = Ai 

. A 4 - 4cos/i Af + (2cos(2/i) + 4)Af - 4cos/i Ai + 1 
Ai = ^ , (27) 

for A = A2 

A 4 . - 4 cos /i A 3 , + (2 cos (2/i) + 4) A 2 , - 4 cos /i A 2 + 1 



Bi 



A 2 -l 



Hence, taking into account Ai • A2 = 1, we have 

B 1 = -^-A 1 . (28) 

Finally, taking account of (26), (27), (28) and |A| = 1, |Ai| < 1, using the 
formula for geometric progression from (25) consequently we get 

F[D]( P ) = 

2 ( ^ 2ftcos(2ft) - sin(2/i) + A + A 1 1 B x 1 \ 



sink — h cos h \ sin /i — h cos /i A ) l-^i A 2 1 — - A - 

\ A A2 

2 ( ^ 2h cos(2/i) - sin(2/i) A 



sink — h cos /i \ sink — h cos /i A 

7=0 ^ ' 7=0 / 
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2 W1 , , , , 2/icos(2fe)-sin(2/i) , A x , 

= £ r • A 1 + Aij H — h - — h 

sin n — n cos /i \ sin n — n cos n, Ai 

OO OO \ OO 

7 — — 2 7—2 / 7— — oo 

Hence keeping in mind A = e 27Ttph we obtain the explicit form (4) of the 
discrete function D[f3]. 
Theorem 1 is proved. 

Theorem 2 is proved using Definition 3 and by direct calculation of the 
left hand sides of l)-5). 

Remark. From (4) we note that D[f3] is the even function, i.e. D[0\ = 
D[—j3] and since \\\\ < 1, then the function D[jS\ is exponentially decreased 
as j3 — y oo. 
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